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Quantum coherence and interference effects in atomic and molecular physics has been extensively 
studied due to intriguing counterintuitive physics and potential important applications. Here we 
present one such application of using quantum coherence to generate and enhance gain in extreme 
ultra-violet(XUV)(@58.4nm in Helium) and infra-red(@794.76nm in Rubidium) regime of electromag- 
netic radiation. We show that using moderate external coherent drive, a substantial enhancement in 
the energy of the lasing pulse can be achieved under optimal conditions. We also discuss the role of 
coherence. The present paper is intended to be pedagogical on this subject of coherence-enhanced 
lasing. 

PACS numbers: 42.65.-k, 42.65.Ky, 42.50.Lc 



I. INTRODUCTION 

Interaction of light with matter is a fundamental ar- 
eas of research in quantum optics and atomic physics. 
Quantum coherence and interference[T| has led to many 
novel effects[2| for e.g. coherent population trap- 
ping ||3]-0, amplification or lasing without popula- 
tion inversion(LWI) [6-9], ultraslow light llOHl2i, en- 
hancement of refractive index without absorption | |T3l - 
[16 1, highly sensitive magnetometry ATI [TSl , coherent 
Raman umklappscattering[19J, high resolution nonlin- 
ear spectroscopypSOl [211, sensing nanoscale molecu- 
lar complexes|22J, photodesorption|23|, bridging quan- 
tum optics with position dependent mass Schrodinger 
equation(PDMSE)|24J etc. Recent experimental and 
theoretical studies have also provided support for the 
hypothesis that even biological systems use quantum 
coherence ll25ti27l ■ Nearly perfect excitation energy 
transfer in photos5mthesis is an excellent example of this. 
Furthermore, during the past decade study of quantum 
interference(QI) effects has been extended to tailored 
semiconductor nanostructures like quantum wells and 
dots due to coherent resonant tunneling owing to their 
potential applications in photo-detection [28, 29 J, las- 
ing |30ll3T]| , quantum computing and quantum circuitry 
1321 l33l , optical modulator|38[. Generally these coher- 
ence are generated using coherent source (laser) to ma- 
nipulate the optical response of the system. But coher- 
ences can also be induced by vacuum also know as vac- 
uum induced coherence (VIC)[34[. Interplay between 
the coherence generated by laser + microwave source has 
also been studied in regard to microwave controlled elec- 
tromagnetically induced transparency [ 35 1, four- wave 
mixing[36J, Raman and SubRaman generation[37J etc. 

Recently quantum coherence effects has extended 
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its domain to plasmonics with a recent proposal of 
coherence-enhanced spaser|39| and propagation of sur- 
face plasmon polaritons[40]. Coherent control in plas- 
monics will definite add a new dimension to the field of 
nanophotonics. On one hand quantum coherence effects 
in quantum optics and atomic physics is a subject of in- 
tense theoretical and experimental investigation while 
on the other hand its effect in human brain has been a 
topic of debate and discussion ll4l1l . 

Coherent excitation in two-level system, studied by 
Mollow[42j, brought interesting features in the reso- 
nance fluorescence spectrum which was later confirmed 
by the beautiful experiment! 43 1. A counterpart of the 
Mollow's triplet|42J was observed with incoherent ex- 
citation in a cavity by Valle and Laussy[44| where they 
showed that the strong-coupling between the cavity and 
the emitter generates the necessary coherence required. 
For multi -level system the coherence can be easily gener- 
ated by coupling the upper-level to an adjacent level with 
a coherent electromagnetic field. Recently Scully ll45l 
extended the idea of coherence effects to solar photo- 
voltaic cells and showed that such devices can benefit 
from quantum boost. In fact this coherence in solar cells 
can be generated by an external source like microwave 
radiation source or by noise-induced quantum interfer- 
ence which is essentially different from the former which 
costs energy ||46ll47l . 

Although numerous theoretical and experimental 
studies of coherence effects have been performed, there 
are still open areas to be explored. For example, quan- 
tum coherence and interference which plays a key role 
in LWI as shown extensively in the literature, the burn- 
ing question we always ask: Can it be used as a tool 
for enhancing the gain in the X-Ray/XUV regimes of elec- 
tromagnetic radiation? A realistic approach in this area 
may open a door for the development of more power- 
ful lasers in the wavelength down to "water window" . 
One approach was proposed by Scully|'48l in which it 
was shown that intense short pulses XUV radiations can 
be produced by cooperative spontaneous emission or 
Dicke superradiance[49| from visible or IR pulses. Later 
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on we also proposed using coherence to generatefSO] and 
enhance [51 1 gain in XUV regime with Helium, Helium- 
like Carbon, Boron as our gain medium. A unique way 
to accomplish effective unidirectional excitation using 
bi-directional source was discussed in Ref.[52J to boost 
gain in the XUV regime. 

In this paper we review and extend the approach of 
applying a strong driving field on an adjacent transition 
to the lasing transition to enhance gain and show that 
gain can be substantially (more than an order of mag- 
nitude) increased under optimal conditions. Here we 
have discussed two regimes, transient and steady-state, 
for coherence enhanced lasing. We have considered las- 
ing on extreme ultra-violet(XUV) transition of Helium 
and Di transition of Rubidium in the transient and the 
steady-state regimes respectively. 

This paper is organized as follows. In section II, we 
discuss the inversion requirement for lasing in the two- 
level emitter based gain medium. In section III we dis- 
cuss the effect of coherent drive origain for three-level 
emitter in A-configuration(see Fig. I21 with initial popu- 
lation inversion(^fl„(0) + QcdO) > Qbbi^)) in the two limits 
of ratio between the spontaneous decay rate on the drive 
transition(yc) and the lasing transition(yi,), (a) yc » yt 
and (b) yc <^ yt- In section IV we consider Rubidium 
laser at Di transition and study the effect of coherent 
drive on the output energy in the steady-state regime. 
In section V, we present the discussion and conclusions. 
We have included appendices to discuss (a) Backward Vs 
Forward gain for three-level emitter based gain medium 
with initial population inversion, (b) A-configuration 
with bi-directional incoherent pump between the lower 
to (dipole forbidden transition), (c) S-configuration with 
uni-directional incoherent pump from lower level \b) to 
uppermost level |c), (d) brief review of vacuum-induced 
coherence between two-levels and (e) brief discussion 
of density matrix vs rate equations for two-level atom 
excited by an external coherent source. 



II. TWO-LEVEL QUANTUM EMITTER BASED GAIN 
MEDIUM 

Let us consider the simplest yet important system in 
laser physics i.e a two-level quantum emitters (semi- 
conductor quantum dots, atoms, molecules, rare-earth 
ions) interacting with a single mode radiation field of 
frequency v. Let \a) and \b) represent the upper and 
the lower levels of the emitter with energy hcOa and Hcoh 
respectively (as shown in Fig.|[Tl). We will study our 
system (emitter+field) in semi-<lassical approximation 
in which we treat the emitter as quantum mechanical 
and the field classically. To describe the the emitter re- 
sponse we will use the density matrix formalism. For 
the two-level medium the evolution of the density ma- 
trix elements Qij takes the form 

Qaa =ra- (^6 + yo)9aa + rgtb - i (oj,Pfl;, - Obgli) , (1) 




FIG. 1: Two-level quantum emitter interacting with a coherent 
field. Field Q;, couples the level \a) and Level \a) decays to 
|b) with a rate jh due to spontaneous emission while there is 
an incoherent unidirectional pumping r from to \a). r„ and 
Tj, are the pumping rate into the levels \a} and \b} respectively. 
yo is the decay rate out of the levels \a) and \b}. 



Qbb = rb+ ybQaa -{r + yo)9bb + i {^iQab - ^bQl^ , (2) 



Qab - -^ab9ab - i^b (Qua - Qbb) , 



(3) 



where yb is the spontaneous decay rate \a) Qb is the 

Rabi frequency for the radiation field. Y^b = yab + ^^b is 
the total relaxation rate of the optical coherence which in- 
cludes spontaneous emission, incoherent pumping, col- 
lisions, detunings etc. Also yab = 70 + {jb + '')/2 + yj^^ 
where y^^^ is the purely phase relaxation and r is the rate 
of incoherent pump from \h) —> \a). fa and rj, are the 
pumping rate into the levels \a) and \b) respectively. }'o 
is the decay rate out of the levels \a) and \b). In the weak 
probe field regime where the population of the levels do 
not depend on the field Q;,. In the this limit, the steady 
state {Qij = 0) populations Qaa and Qbb are given by 



rn{r + yo) + nr 

yo{r + yb + yo)' 



(0) _ niyb + yo) + rgyb 
" yo{r + yb + yo) 



Solving for Qab we obtain. 



„ _ (0) on 

Qab --^TT- [Qaa ~ Qbb } ' 
- _ nC) 



(4) 



(5) 



(6) 



The inversion defined as W = gj - q\^ is not reached if 



Q^'m < Q^l which requires 



yo{ra - n) < {yb - r){ra + n). 



(7) 



When we do not consider any pumping into or out of the 
level, i.e Va, Vb, yo = 0, and the population is conserved 
i.e Qaa + Qbb = 1/ the emitter does not show population 
inversion in the steady-state if incoherent pump rate r 
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is less than the spontaneous decay rate yi,. 
definition of linear susceptibility we obtain 

(1) . 3 3 \(^aa fjf,; 



Using the 



(8) 



The imaginary part of the complex susceptibility is given 
by 



(9) 



The weak field Q;, will be amplified if Im;^' < 0. Thus 
from Eq.||9}, we obtain the necessary condition for gain 
as 



kfaa ^ tfbb ' 



(10) 

Equation JTOl is also known as population inversion 
condition for two-level system. It is worth mentioning 
here that Mollow garn|53| (or hyper-Raman, or three- 
photon gain) can be obtained for two level system even 
in the absence of population inversion in the bare ba- 
sis. In fact this process does not require any popula- 
tion in the upper level. Here two pump photons are 
absorbed while a probe photon is emitted. Mollow 
gain is a Raman-like process where the energy is trans- 
ferred from the pump to the probe beam and thus it 
is different from amplification without inversion (AWI) 
where we are interested in extraction of energy from 
the medium. For phase effects in these hyper-Raman 
process in effective two-level system and using phase 
jump to control the excitation the readers are suggested 
the experimental] 54, 55 1 and theoretical|56, 57 1 papers 
respectively. It is worth mentioning here that exact solu- 
tions for the transient probability amplitudes Cafi{t) for 
two-level atoms interacting with ultra-short pulses can 
be obtained analytically Il58ll59l . 



III. COHERENCE-ENHANCED LASING I: TRANSIENT 
REGIME 

A. Three-level quantum emitter based gain medium 

In the previous section we briefly reviewed the pop- 
ulation inversion condition required for lasing in two- 
level emitter based gain medium. If we add another level 
to the emitter, the physics of the light-matter interaction 
becomes rich and has opened the door for many coim- 
terrntuitive physics with vast application. In this section 
we will review the concept of coherence-enhanced lasing 
in transient regime as discussed in Ref.[51 1. 

Let us consider the three-level quantum emitter in 
Lambda (A) configuration. Here for the sake of sim- 
plicity consider the gain medium as three-level atoms in 
which the transitions \a) |c) and \a) \b) are electric- 
dipole allowed but the transition |c) <r-> \b) is electric- 
dipole forbidden [see Fig. [2| due to selection rule based 




FIG. 2: Three-level atomic system in (a) A and (b) H configura- 
tions. 



on parity. We will assume that at the initial moment of 
time the population is distributed between levels \a) and 
It)) onlyi.e p„fl(0) + ^)bi,(0) = 1. Any initial population in the 
level |c) is not a part of the gain medium in the absence of 
the drive field. Transition \a) |c) is driven in resonance 
with the Rabi frequency Qc. We investigate how a weak 
laser seed pulse at the \a) \b) transition evolves during 
its propagation through the medium. The Hamiltonian 
in the interaction picture can be written as 

^ = Ac\c){c\ - [(Qb \a){b\ + Q, |fl>(c|) + H.c], (11) 

The spontaneous decay in the channels ac and ab are 
quantified by the rate yc and y;, respectively. Incorpo- 
rating these decay rates, the equation of motion for the 
atomic density matrix is given as 



(12) 



where, 

Ob = \b){a\ , al = \a){b\ , and a, = |c)<fl| , = |fl)<c| , (13) 

Here Q/, and Qc are the Rabi frequencies of the probe and 
drive fields respectively. Evolution of the atomic density 
matrix gij is described by the set of coupled equations 121 

gab ~ -^abQab - i^bigm - gbb) + i^cgcb, (14) 
gcb = -Tcbgcb + K^'cQab - ^bQ'ac)' (15) 
gac = -Tacgac " i0.c{gaa - Qcc) + i^bQ'cb' (16) 

gbb = ybgaa + i {^Igab - ^bgl^) > (17) 



gcc = ycgaa + K^'cQac - ^cglc)' 



(18) 
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Qaa + Qhh + Qcc = 1, 



(19) 



Qtth{i, z) 



iAi,t-iki,z 



(30) 



where = {y^ + y/,)/2, = {yc + yb)/2 + iAc, Tcb = -/A,, 
are the relaxation rates of the off-diagonal elements of the 
atomic density matrix. Here we assume that the spon- 
taneous decay rate on the drive transition {\a) |c)) is 
much larger than the probe transition {\a) \b)). In this 
limit, the coherent drive will redistribute the population 
between the two-levels \a) and |c) and induce a non-zero 
coherence between them. For simplicity we will assume 
that the drive field is so strong that any variation can 
neglected i.e =const (a real number), Q;, is very small 
and the drive transition is excited resonantly (Aa = 0). In 
this limit we can find analytical expression for the coher- 
ence g„c and the populations gu. The equation of motion 
is given as (putting yy = 0) 



ycQaa + i0.c{gac - C.c), 



Qac ~ ^ acgac ^^cigaa gcc)- 



yaa ' fee 



Pafl(O), 



(20) 



(21) 



(22) 



The steady state solution (gij) can be easily obtained 
asHSll 



4Q2 



yi + 8Qi 



(23) 



(24) 



(25) 



where gaa{0) is the initial population of the level \a). Evo- 
lution of the weak laser pulse Q;, is described by the 
coupled Maxwell-Schrodinger equations 



gab - --^abgab " i^bigaa - gbb) + i^cgcb (27) 



gcb = i{0.cgab - 0.bglc)- (28) 

where rjaj, = {3/8n)NA^j^yi, is the coupling constant, N 
is the atomic density and Aab is the wavelength of the 
\a) <-> \b) transition. Here we also assumed that the drive 
transition has been excited resonantly (A^ - 0). Next 
we move to find the dispersion relation by looking for 
solution of Eqs. | 26p8| in the form 



Qfc(f,z) ~ e' 



fAijt-ikiz 



(29) 



gcb{t,z) 



(31) 



which yields 



iyAb 



. {ckt - Ah) + cAbijabigbb 
^ / (32) 



+CTJab^cgac = 0, 

here A;, is the detuning of the laser pulse frequency from 
the \a) \b) transition frequency. Now the imaginary 
part of kb determines the gain (absorption) and takes the 
form 



lm{kb) = i]ab 



ycAligaa - gbb)/2 + Oc (n? - Afj lm{g„,) 
(A2-Q?)' + (y,A„/2)2 



(33) 

In particular, for the mode resonant with the \a) 1^7) 
transition A;, = and we obtain 



G = lm{kb) = ^lm(pac)- 



(34) 



When we look at the expression for g„c, according to Eq. 
1 25 1, lm{g„c) > which implies that we have gain in this 
configuration even when we have very little population 
in the upper level \a). Next we will study the regime in 
which the spontaneous decay rate on the drive transition 
(|a) |c)) is much smaller than the probe transition 
(|fl) —> \b)). In this limit we have to consider the transient 
behavior of the populations and coherences to study the 
gain. 

Before we discuss this scenario let us briefly review 
the essence of gain in transient regime. The equations 
of motion for the density matrix element g„b is given by 
Eq.|[T4||, which in the steady-state gives 



^ab^lgab] = -iOtigaa - gbb) + i^^cgcb 



(35) 



where gjj is the steady-state value of the density matrix 
elements gij. Now let us draw some general conclu- 
sions from the density matrix equations. From Eq.dSSl it 
appears that, for two-level model i.e Qc = 0, ampliHca- 
tion condition (3 [pa;,] < 0) requires population inversion 
gaa > gbb- On the other hand for three-level system for 
sufficiently negative Qc gcb, the necessary condition for 
the probe transition \a) <-> \b) to exhibit amplification can 
be satisfied even without population inversion. If we 
look at the evolution equation for the population in level 
\b}, 



2Qb^lg^ 



>abi 



ybga, 



dgbb 
dt 



(36) 



Amplification condition (3[^)flb] < 0), reduces Eq.|36 ' to 

Ybgaa < (37) 
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0,,(O) = 0.55, e„„(0) = 0.45 
0„(O) = 0.50, e„„(0) = 0.50 
0„(O) =0.45, e„a(0) = 0.55 
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FIG. 3: Plot of the gain G„(,(t) as function of normalized time jct. 
Solid line (Blue) is for initial Raman Inversion i.e Qcc{Q) > QtbiO) 
while the dashed line (Green) is for zero initial Raman inversion 
i.e [JbbiP) = QcciO) and dashed dot line (Brown) is for initial no 
Raman inversion i.e QcdO) < QbbiO). The common parameters 
taken are fc = l,Yb = 0.2, = 3,Qf, = 0.01,g„„(0) = 0. For 
initial Raman inverted medium QcdO) = 0.55, gbb{0) = 0.45 and 
Qa-{0) = 0.45, ^)(,(,(0) = 0.55 for non-inverted case. We have 
normalized the rates and Rabi frequencies with respect to the 
decay rate on the drive transition \a) — > |c). 



From Eq.|[37| we conclude that in transient regime, the 
probe transition exhibits amplification when the tran- 
sient growth of the population in level \b) exceeds due 
to the incoherent radiative decay \a) \b). Thus a net 
amplification can be realized in during the time of inter- 
action t iff L60J 



Qbbit) > QbbiO) - Yb Qaa{t')dt' 



Jo 



(38) 



In steady -state Eq.|[37[| can never be satisfied but it can 
be realized in the transient regime I(50ll60l . To quantify 
the gain(absorption) profile for the lasing transition, we 
analyze the parameter Gab{t) defined as 



8n 



(39) 



Here N is the number density of atom (ions) and A;, 
is the wavelength corresponding to \b) transi- 

tion. In general Gab{t) is an oscillatory function of t. 
When Gab{t) > 0, the probe pulse will enhance while for 
Gab{t) < it will experience attenuation. In Fig. |3j we 
have plotted Gab{t) for three choices of Raman inversion 
(positive, zero, and negative) respectively. 

In the limit y;, » we can find approximate analytical 
expression for the temporal evolution of the coherence 
Qac and the populations Qij for a weak probe field Q(, and 
strong drive field. We obtain the solutions as 



gaa{0) {sin^(QcO + cos(2QcO 

_yb_ 

4Q, 



sm{2Qct)} , 



(40) 



O 

> 



0.5 
0.0 
-0.5 
-1.0 



•i / \ 

\1 ^' 







1.1 



2.2 



3.3 



4.4 



Time, t (ns) 



FIG. 4: Inversion W(t) in the probe transition (a <-> h) vs time t. 
Dashed black curve shows the result for Qr = 5vb while solid 
red curve is obtained with no drive. For numerical simulation 
we used = 1-1 x lO'^yj and the initial condition p„„(0) = 0.8, 
QaiQ) = 0.2, p,,(0) = g„M = g„t(0) = 0. 



sin(Ocf) 



(41) 



(42) 



■ cos(QcO} ■ 



Defining the population inversion W{t) = Qaa{t) - Qbb{t) 
we obtain the expression 



W{t) = 



QaaiO) 



-yht/2 



[3 + cos(2QcO- 



Or 



sin(2Qcf) 



(43) 



1. 



In Fig. |4] we plot the population difference W(i) as a 
function of time for initial conditions ^)„fl(0) = 0.8, QbbiO) - 
0.2 and Qcc{0) = 0. Solid line is obtained for - 5yb 
while for dashed line = 0. Driving the \a) |c) 
transition yields oscillations in the population difference 
between \a) and \b) levels. The oscillations are due to the 
population oscillating between the levels \a) and |c) in 
the presence of the coherent drive. 



B. Helium as the gain medium 

Next we solve Eqs. |[T4)-|[T9) and p6| l numerically and 
obtain evolution of the proEe laser pulse 0;,(f,z) when 
the \a) <r^ \c) transition is driven by a constant coherent 
field . We perform simulations for the initial condition 
QUO) = 0.8, QbbiO) = 0.2, p,,(0) = and take i]/yb = 40.75 
cm~^ and = 1.1 x 10~''yb. As an example, we consider 
He for which states 2^So (|c)- level), 2^Fi (|fl)-level) and 
the ground state I^Sq (|f')-level) form A-scheme [see 
Fig. p). The model parameters are A„i, = 58.4 nm, Aac = 
2059 nm, y<: = 2 x 10^ s'^ and yb = 1.82 x 10^ s'^ Then for 
ion density N = 10^® cm"'' we obtain rj/yb = 40.75 cm~^. 
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Yc=2x 



Yk= 1.82 X loV 




^c = 728.3 nm 

1.83x lo's' 

a — 



.3's 



\= 1.82 X loV' 



Lasing transition - 



-2'P, 



58.4 nm 



-I'Sn 



FIG. 5: Energy level diagram of He atom in A (left) and H 
configurations (right) respectively. Lasing at 58.4nm can be 
enhanced using a coherent drive at 2059nm(A-configuration) 
or 728.3nm(5-configuration) l5Tl 



We assume that input probe laser pulse has a Gaussian 
shape 




2 4 6 8 
Propagation length, L (cm) 



FIG. 6: Ratio of the output energy to the input energy of the 
probe laser pulse as a function of sample length L with three 
choices of drive Rabi frequency Qc = 0,9 and 12.6 x lO's"^ In 
numerical simulations we take fc = 1.1 x lO'^yt, Tilyt = 40.75 
cm~^ and assume Gaussian initial probe pulse shape given by 
Eq. |44|. Initial populations are Qaa{G) = 0.8, QbbiP) = 0.2 and 
Qcc(0) = 0, while initial coherences are equal to zero. 



Cli{t,z = 0)= 0.01 exp 



ytt- 0.12 
0.05 



(44) 



During propagation of the weak laser pulse through the 
medium the atomic population spontaneously decays 
into the ground state. After a certain time the medium is 
no longer inverted and the laser pulse begins to attenu- 
ate. Thus, there is an optimum length of the atomic sam- 
ple which yields maximum enhancement of the pulse 
energy. For the optimum length the pulse leaves the 
medium at the onset of absorption. In Fig.|6]we plot the 
ratio of the output pulse energy to the input energy as a 
function of the sample length for three choices of exter- 
nal drive. We find that optimum length corresponding to 
maximum output energy without any drive is Lq - 4.25 
cm. This optimum value depends on the initial popu- 
lation inversion and the decay rates. At this optimum 
length the ratio of the output to the input probe field 
energy is ~ 1.42 X 10''. In the presence of drive field 
Qc ~ 12.6 X 10's~^, this ratio at the propagation length 
(Lq) optimum for two level configuration increases to 
5.81 X 10^. Hence we see that coherent drive can increase 
the energy of the XUV lasing field by 4-fold. To optimize 
this enhancement with respect to the drive field, but keep 
the sample length to be Lq - 4.25 cm, we have plotted 
the ratio of the output to input energy as a function of 
Rabi frequency of the drive field. However its worth to 
mention that this length Lq does not corresponds to the 
maximum gain for three-level system as clear from dot- 
ted circles which corresponds to optimum propagation 
length in Fig.|7] 

In Fig. [7] we plot the ratio of the output laser pulse 
energy (at z = Lq) to the input energy (at z = 0) as a func- 
tion of strength of the driving field Clc- One can see that 
in the presence of coherent drive the output pulse en- 
ergy oscillates as a function of Oc and for some optimum 



value of Clc we observed 4-fold enhancement with drive 
with respect to no drive. The coherent drive change the 
dynamics of the system in two ways (a) it redestributes 
the population between the upper two levels and (b) 
induce Raman coherence Qcb which is the key to any las- 
ing without inversion(LWI) schemes. For < Qc < 2, 
the drive effectively acts as an incoherent pump i.e de- 
populating the upper level \a) \c) and the ratio goes 
down from 1.41 X 10^{Clc = 0) to 4.84 x 10^^ = 2). 
After Qj. > 2 the ratio starts to increase and goes back 
to the value of 1.42 x 10^ at - 4.625. Beyond this 
value of Qc, the ratio further increases and reach the 



Y,,ft=0 V=1.8 Yp,,=3.6 YpA=5.4 YpA=7.2(xl0^s-l) 
6 F 




FIG. 7: (a) Ratio of the output energy to the input energy of 
the probe laser pulse as a function of the driving field Rabi 
frequency Qc. The ratio is ~ 1.42 x 10* at Qc = 0. In nu- 



1.1 x 10-3}'6, ri/yt = 40.75 



merical simulations we take y, 
cm 

Eq. _ 

the initial populations are g„a{0) = 0.8, gbb(0) = 0.2, Pcc(O) = 
and pac(O) = QabiO) = 0. Other curves corresponds to different 
choices of the dephasing rates = 1.8,3.6,5.4and7.2xl0**s*^ 



^ and assume Gaussian initial probe pulse shape given by 
JiiJ. The length of the sample is L = 4.25 cm, while 
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global maxima 5.80 x 10'' at Qc ~ 7.25 thus an enhance- 
ment of 4-fold with respect to no drive Qc = 0. For 
7.25 < < 14.625 the ratio shows similar behavior with 
the local minima of 4.19 X lO^Qc = 10.625) and local 
maxima 3.78 X 103(Oc = 14.625). Thus, coherent drive 
can increase the laser pulse output energy as compared 
to the pulse energy with no drive. We have also plotted 
the ratio when we phenomenologically add a dephas- 
ing term (due to collision, rnhomogeneous broadening 
etc) to the coherence for different values of ypu ■ We see 
that the enhancement factor decreases with increasing 
dephasing rate and at ypf, - 10's~^ we do not see any 
enhancement due to drive. Thus, one has to cleverly 
reduce the dephasing rate to observe the nice coherence 
enhanced lasuig in XUV regime. 



Qaa + Qbb + Qcc - 1/ 



(53) 



where Tab = (Yc + yb)/2,Tac = {fc + }'b)/2 + iAaJcb = 
y„/2 + iA„ are the relaxation rates of the off-diagonal el- 
ements of the atomic density matrix. The steady state 
solutions of the off-diagonal elements of the density ma- 
trix is calculated from Eqs. | 4&pO| by setting the time 
derivatives terms to zero. Solving we get. 



. . j n„i,{\D.b\'^ +TcbTca) + nJQaf- 



Qch = 0.tCla 



ca^ ah 



Pfcpr^b + (TcfcT^b + |Q«|2)r„ 



(54) 



(55) 



IV. COHERENCE-ENHANCED EASING II: 
STEADY-STATE REGIME 

In this second part of the paper, we will study 
coherence-enhanced lasing in steady-state regime. We 
will begin with the generic three-level scheme as shown 
in Fig. l2|b) i.e emitter in the Cascade/Ladder configura- 
tion. Let us also assume an incoherent pump from the 
ground level \h) to the intermediate level \a). Transition 
\c) <r-> \a) is driven in off-resonance by the drive field. The 
Hamiltonian in the interaction picture can be written as 

^ = Aa\c){c\ - [(Q„ \a){b\ + \a){c\) + H.c], (45) 

The spontaneous decay in the channels ac and ab are 
quantified by the rate ya and yb respectively. Incorpo- 
rating these decay rates, the equation of motion for the 
atomic density matrix is given as 



g = -iYV, + y ([Ob, pal] + [Obp, a\;[) 



(46) 



where, 

Ob = \b){a\ , a\ = \am , and ct„ = |fl)(c| , a\ = \cW\ , (47) 

Here Q;, and □« are the Rabi frequencies of the probe and 
drive fields respectively. Evolution of the atomic density 
matrix is described by the set of coupled equations fJ] 

Qab = -TabQab - i^biOaa - Qbb) + i^'aQcb, (48) 
&cb = -'CcbQcb + i^aQab - i^bQca, (49) 
Qca - -^caQca " iCiaiQcc " Qaa) " i^lQcb, (50) 



Qbb = ybQaa + I 



(plQnb - ObQlb) , (51) 



WflfcPbp + UcaiTcbTab + PP) 

Pbpr„b + (r,6r„6 + p„|2)r,„ 



(56) 



Using the condition Qf, « Oa, further reduces Eqs. 1 54 
|56lto. 



'n^XbTca+nfJlQf' 

{T,bTab + P«P)r.« 



{T,bT,b + P«P)r„ 



(57) 



(58) 




2 4 6 
Detuning, A a 
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Qcc = -yaQcc - im,Qca - ^aQla), 



(52) 



FIG. 8: Plot of -Im[^„i,]/Q[, as a function of the detuning A„ 
for different values of the drive Rabi frequency In (a) we 
have we assume incoherent pump from the lower level I fe) —> |c) 
while in (b) we considered the incoherent pump from the lower 
level —> \a}. The dashed line corresponds to no drive Q„ = 0. 
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a^i' = -/Q„|^| (59) 

where - g^'' - Here, g^^'' is the population of the 
level |f) where we keep the probe field Qj, to its lowest 
order while we keep all the terms for the drive field 
Here we have considered to scenario of unidirectional 
pump from (a) \b) \c) and (b) \b) — > \a). 

In Fig.ljSJ we plot the imaginary part of gab for the 
two scenarios mentioned above. When we incoher- 
ently pump to the upper level |c), we observed gain 
enhancement when we apply moderate but not strong 
drive. For eg when = 1, and = 0, the gain oc 
-lmp„;,/Ofc is enhanced by 6-fold. By introducing detun- 
ing on the drive transition this enhancement decreases. 
When 15/8 < Aa < 33/10 the gain switches its sign and 
we observe loss. Thus the gain medium becomes absorp- 
tive for appropriate range of detuning Aa and can be tune 
back to show gain beyond this range. For A^ > 33/10, 
transition \a) <r-> \b) shows gain and reaches the value of 
gain with no drive at large detuning. Similar behavior 
(qualitatively) is also seen at other drive Rabi-frequency. 
On the other hand when we apply incoherent pump 
to the level \a) and the drive field is detuned, the gain 
with Qfl is lower than no drive as shown in Fig|8|b). 
When the detuning is further increased the gain reaches 
asymptotically to no drive gain. The dashed line in Figs. 
|8] corresponds to gain without drive, which should be 
immune to detuning as shown by constant value while 
varying A^. 




FIG. 9: (a) Energy level diagram of atomic Rubidium. Here we 
have the dipole allowed transitions 5^Si/2 <-> 5^Pi/2 (Di)line 
and 5^Si/2 <-> 5^P3/2 (D2)line. The population between the 
levels Pi/2 and P3/2 are exchanged due to collisions by buffer- 
ing the alkali vapor with other gasses like helium, ethane etc. 
(b) Four-level model for coherence enhanced Rubidium laser. 
Here couple the drive transition with a coherent field of Rabi 
frequency The Rabi-frequency of the pump and lasing 
fields are CI,, and Q;, respectively. 

via collisions with the buffer gas. In this paper we 
have considered ethane as the buffer gas for which we 
evaluated the collision rate as R„e = 4.383 x 10^s~^ and 
Rea = 3.245 X lO^s"^ After ~ 10ns the system reaches 
steady-state with population inversion on the lasrng 
transition as shown in Fig.|[TO|. The steady -state pop- 
ulation inversion is fiab - 0.06but the overall population 
inversion i.e psn + gee - gtb - 0.369. 



A. Rubidium Laser 

The level structure for the Rubidium laser (Djline) is 
shown in Fig. l9] In the last section we discussed coher- 
ence enhancedlasing in Helium in the transient regime, 
here we will show that in the presence of the drive field 
Ofl gain on the lasing transition (Dilrne) can be enhanced 
substantially even in the presence of a strong dephasrng 
rate ypi,. Detailed analysis of Rubidium laser, and the 
conditions under which lasrng action can be achieved 
on (Di) transition ha been discussed in|61 1. 

One important condition to achieve population inver- 
sion on the lasing transition is that the rate of popula- 
tion exchange between 5^^3/2 and 5^Pi/2 should be much 
faster than the rate of spontaneous decay from the level 
5P 5S. The transition 5^P3/2 5^Pi/2 is electric dipole 
forbidden, hence the population exchange is achieved 
using collisions with buffer gas. Here in this section we 
have assumed He as the buffer gas. The excitation trans- 
fer cross-section for Rb induced by collisions with rare 
gas atoms and alkali metals can be found in[62J. 

We have modeled our gain medium as a four-level 
system as shown in Fig. Mb). We excite the transi- 
tion |e) \b) hy a pump laser and the population 
between the dipole forbidden transition is exchanged 



B. Steady State Gain 

The equation of motion for the density matrix elements 
gij are given as, 

gab = -Tabgab " i^bigaa - gbb) + i^lgcb - i^eg'ea' (60) 

gca = -^cagca - H^aigcc - gaa) - i^'bgcb, (61) 

Qeb = -'^ebQeb - i^ei^ee ~ Qbb) - i^bQea, (62) 

gcb = -^cbgcb + i^agab - i^bgca - i0.egce, (63) 

gee = -Tcegce + i^agla " i^'eQcb, (64) 

Qea = -TeaQea + i^eQ^b " ^^bQeb " i^aQ'ce- (65) 

The population terms is given as 

Qaa = - (Yb + Rae)gaa + Reagee + Yagcc ^^^^ 
+ K^'aQea - ^aQ'ea) " K^'bQab " ^bQ'^b)' 
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FIG. 10: Temporal evolution of the populations of the level \a) 
(solid black line) and lower level lb) (dashed red) of the lasing 
transition in the absence of the drive field Q„. For numeri- 
cal simulation we used Rga = = 0.74:, ft = 0.085, = 
0.087, y„ = 0.0085, Q(, = 0.001, Q, = 5. Here we have normal- 
ized the decay rates and the Rabi frequencies to the population 
transfer rate R„, which is an order of magnitude higher than 
the spontaneous decay rates. 



Pbb = ybQaa+}'eQee + i{0.lQab-0.bQli^) + i{0,lgeb-Oegl^), (67) 



Qcc = -JaQcc - im,Qca - ^aPla), 



(68) 



FIG. 11: Plot of steady-state gain G„i, as a function of the co- 
herent drive Q„ for different values of the dephasing rate Vph- 
Other curves corresponds to different choices of the dephas- 
ing rates fph = 0.0,0.5,1,2,10. Other parameters are same as 
;IT0l 



(75) 



It can be easily verified that we can obtain the known 
results p , for cascade and Vee scheme. The zeroth order 

^ab . II . 

population obtained from Eqs. 1 66|69| as 



Qee - -{ye + ^ea)Qee + KeQaa - K^lQeb - ^eQ^^)/ (69) 

where. 



r 

J- ca 
T 

^ ea 



Ya + yb+ Rge _ 7 a + Ye + 
^ / ^ cc ~ Z 



ye + yb+ Rge + Rm y _ 7^ 
2 ' ^ cb — 

Ye + Rm _ yb+ Rea 

/ Iflfc- 



(70) 



(0) _ 2Rea{IaIea + 2Q„^)Qg 
Qaa 



(0) ^ 2R,„QgQf 



(76) 



(0) plReaYbTeb + (Roe + Yb){ye^cb + 2Q?)] 

ib = }vi ^^^^ 



(78) 



Let us assume that all the fields are real and we keep the 

probe field Q;, to the lowest order while we consider all 

orders for the drive fields and [see Fig. l8|b)]. In 

(■\\ ' — ' 
this limit, we obtain the coherence g^^ as 



„(0)^ + „(0)S_„(0)C 

ab '■" eb 



C - Tea \TcbTce + Op - Q^j Op 



(71) 



where, the population inversion term = - pj.^' and 
the constants are defined as 

= TcflTrf, [rcb(rcerpa + O^) + TeaOp j 

s = r,b [r,j,„ + o2 - o2] o2 



(72) 



(73) 



(74) 



2{Ra, + yb)0^ 
M 



(79) 



where the constants /3 = {y„Tc„ + 20^) and M = 
YaTca \Rea{ybTeb + 20^) + (R„, + yfc)(y,r,fc + 40^)]. In the 
absence of the drive field Oa, and strong pump field 
Of » Ye, Yb,Rm,Rae we obtain 



Re, 



Rea+1{Rae + Yb)' 



(0) _ Rae + Yb 

Qbb 



Rea+2{Rae + Yb)' 



JO) _ 



Rge + Yb 
Re,+2(Rae + Yb)' 



(80) 



(81) 



(82) 
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FIG. 12: Plot of steady-state gain G„;, as a function of the coher- 
ent drive Q„ for different choices of the detuning A„ = 1, 2, 3, 4. 
Here A„ normalized to the population transfer rate R^a- Other 
parameters are same as Fig. [T0| 



The steady-state inversion {gfj + gf^ - g^^') is given as. 



Rea + 2(Rae + Yb) 



>0 



(83) 



and also on the lasing transition {a <-> h) we obtain. 



a. 



Rea + 2{Rae + YbY 



(84) 



For inversion on the lasing transition we require Rea > 
Rae + Yb- This is a very crucial condition to keep in mind 
while selecting appropriate buffer gas for population ex- 
change between the upper levels |e) and \a). We have 
also included here an appendix to review similar con- 
ditions for three-level emitter in Lambda-configuration 
with symmetric bi-directional pump for the convenience 
of the readers. 

Let us now study the effect of drive field Cl^ on the 
steady-state gain Gab defined b y Eq.p9| . The result of the 
numerical simulation of Eqs. | 60||69| is shown in Fig. 11 
in which we have shown the eHect of the drive field on 
the gain Gab- We see that in the presence of the drive 
field we can enhance the gain by an order of magnitude 
for Cla ~ 5. To emphasize the role of the coherence 
we simulated Eqs. | [60p9) in the presence of the several 
dephasing rates fpi, andthe results are shown in Fig 1 1 
We see that even when the dephasing rate is 20-foIc 
higher than the spontaneous decay rate in the lasing 
transition, quantum coherence can still enhance the gain 
by 2-fold. 

In Figll2]we have plotted the steady-state gain Gab de- 
fined by Eq. 1 39 1 as a function of the coherent drive Qa 
for different choices of the detuning A„ = 0,1,2,3,4. In- 
terestingly the central peak splits into two peaks which 
further move apart as we increase the detuning. The first 
peak vanishes beyond = 5. The position of the peak 
is also asymmetric with respect to the peak with zero de- 
tuning (Aa - 0). We also observed that within the range 
of parameters (drive detuning and drive Rabi-frequency) 



used to simulate Fig. (12) there was no significant change 
in the width of the second peak. 



V. CONCLUSION 

In this paper, we have studied one of the manifesta- 
tion of quantum coherence effects in atomic physics. We 
showed that quantum coherence and interference effects 
can be used as a tool to enhance gain in both XUV and 
infra-red wavelength under suitable conditions in tran- 
sient and steady-state regimes respectively. We investi- 
gated the role of coherence in substantial enhancement, 
more than an order of magnitude, in the output energy 
when compared to no drive model (Qa = 0). 

For transient regime we selected Helium atoms (lasing 
at 58.4nm), initially prepared in the highly inverted state, 
our gain medium. We optimized the propagation length 
for a given population inversion by studying the energy 
of the seed pulse with varying propagation length. At 
the optimum length, the ratio of output to input en- 
ergy is 1.41 X 10^. At this optimum length we applied 
a coherent drive (Aa ~ 2/.im ) and optimized the drive 
Rabi-frequency to enhance the gain. For moderate drive 
we theoretically demonstrated a 4-fold enhancement in 
the ratio. For steady-state regime, we selected Rubid- 
ium atoms (lasing at 794.76nm) as our gain medium. 
Here we observed more than an order of magnitude en- 
hancement in the steady-state gain even when the phase 
relaxation rates (dephasing) along with other decoher- 
ence rates are 20-fold faster than the spontaneous decay 
rate on the lasing transition. These results though seems 
optimistic in the XUV regime, which is the main moti- 
vation for this project, can provide a route to building 
powerful laser if we can cleverly handle and reduce the 
fast relaxation rates. A proof of principle experiment on 
coherence-enhanced lasing with population inversion using 
favorite atoms like alkali metals (Rb, Na, K) would be one way 
to go. 

With the recent proposals[39 40] on extending the co- 
herence effects to new domain of plasmonics has defini- 
tively opened new prospects and challenges which has 
been unexplored till now. Moreover optical control or 
rather coherent control in plasmonics will add a new 
dimension in nanophotonics. 
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Appendix A: Backward Vs Forward Gain 

In this paper have considered the evolution of the in- 
jected seed pulse at z = in the forward direction (with 
respect to the driving laser) . In this section we will briefly 
discuss the evolution of an identical seed pulse injected 
at in the backward direction at the end of the sample i.e 
z = L along with the forward seed pulse. We will con- 
sider the three-level system in Lambda configuration in 
the limit yi, » as shown in Fig. l2] and we drive the 
transition a c by a field E^. We also excite the transi- 
tion a bhya weak probe field Ei,. We write the electric 
field as, 

£,(z,0 = f [£,"(z,fy^'^+c.c], (Al) 



E^(z,0 = ^[£,^(z,f)e'^?+c.c], (A2) 



Using the transformation Eqs. | [A7|A10[ | in Eqs. | [A4]|A6 i 
we obtain for the backward direction: 



where 



- kz- vt, 



= -kz - vt 



(A3) 



Here (+) and (-) sign as the superscript means forward 
and backward direction respectively. We can write the 
off-diagonal term as 

pah = -COabpab - ipab ' Efc(Pflfl - Pbb) + ipac ' Ecpcb (A4) 



Pac = -(^acpac - ipac ' E^paa - Pec) + Wab ' Etp*^ (A5) 



Pcb - -(^cbpcb - ipab ' EfePac + Wca ' Ecpab 

Let us make the following transformation 



is* 



Pac = QacS' 



Pcb = 9^"*' + q;/""'' 



where 



(A6) 



(A7) 



(A8) 



(A9) 



= e--dt (AlO) 



Qab = ~^ab9„b " ^^b (O"" " Qbb) + i^cQ.i, 
Q~cb = -^cbQ'^ - i^bQac + i^'cQab 

and for the forward direction we obtain 

9ab = -^"bQ'^b ~ i^t(9an " Qbb) + i^cQ'^b 
Qac = -'^acQtc - i^ciQaa - Qcc) + t^t Q^h 

The evolution of the population is given as 

9aa = - {yb + yc)eaa - W 9ab " "fc P«fc) " ^i^b' &ab 



(All) 
(A12) 

(A13) 
(A14) 
(A15) 

(A16) 



gtb = ytoaa + iicifg^, - n^g^;) + f(Q-*^)- - q-^^-*) (ai7) 



Qcc = JcQaa + K^TqIc ' ^UVc) 



(A18) 



From Eqs. | [AT6]|AT8) we see that the population equa- 
tions are symmetric under the transformation + -, 
hence the evolution of the injected (identical) seed pulse 
at the respective ends of the sample will be the same. 

If we can selectively destroy the coherence in the back- 
ward direction (^Jj^), keeping the forward coherence in- 
tact, we can in principle get asymmetric results as dis- 
cussed in [51]. In the opposite limit, than the results 
reported in [51 J, greater energy in the backward direc- 
tion will be beneficial in the sky laser physics [63]. 



Appendix B: Three-level lambda with bi-directional pump 

Let us consider three-level emitter in Lambda (A) con- 
figuration as shown in Fig. 13. Here we have included 
a bi-directional pump given by the rate r to exchange 
population between the lower two levels. The evolution 
of the density matrix elements is described by the set 
of coupled equations Q: 

gab ~ -^abQab + i^bi^bb - Qaa) + i^cgcb, (Bl) 



Qcb = -'^cbQcb + K^'cQab - ^bQac)' (^2) 



gac = -^acgac - i^ciOaa - Qcc) + i^bPcb' (^•^) 



Qbb = rQcc - rgbb + YbQaa + i(nigab + n,gl,), (B4) 
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FIG. 13: Three-level atomic system in A— configuration. Here 
we have ir\cluded a bi-directional pump given by the rate r to 
exchange population between the lower two levels 



Substituting Eqs.| [B7p9| l in Eq. | |B12[ | gives. 



(r,i,r„fc + Q2)Q 



(B13) 

where k = -{yt + Yclf^ca- In the strong drive field limit, 
to observe gain m the probe transition requires Im^af, < 
requires 



Yc > Yh, 



(B14) 



which means the rate of spontaneous decay in the driv- 
ing transition \a) \c) should be greater than probe 
transition \a) \h). Thus to observe gain without pop- 
ulation inversion the rates should satisfy Eqs.( Bll) and 
I B14| simultaneously. 



Qcc = -rQcc + rQbb + YcQaa + K^'cQac - ^cQIc)> (B5) 



Qaa + Qbb + Qcc = 1- 



(B6) 



where r,y = r..,Yab = Tac = (Yb + Yc)!'^ and Ycb = r. Here 
Qfc and Oc are the Rabi frequencies of the probe and 
driving field respectively. r,y is the relaxation rates of the 
off-diagonal terms Qij. Similar to the two-level system, 
we obtained the steady-state (long time limit) solution 
for the populations Qu in the weak probe field and Oc is 
real. 



2Q|r 
Q ' 



(0) (Yb + Yc)r^ca + 2{yb + r)D^ 
&bb = 



- 



Q 



(Yb + Yc)rTci, + 2?-Qg 
Q 



(B7) 



(B8) 



(B9) 



where Q = 2r(Yb + Yc)^ca + 2{Yb + 3j')Qc. Population 
inversion can be calculated from Eq.(18-20) as. 



2{r-Yb)^l 
Q 



(BIO) 



FromEq.|[BTO) we obtain the no population inversion 
165] [gfj +^^- gfij < o) condition as 



Yb > r. 



We obtain gab as 



gah-'^b' 



[fan fj,;, J ^ cb^ ac ^ [fee fua J '"'i: 



{^cb^ab + ^a, 



(Bll) 



., (B12) 



Appendix C: Three-Level E-configuration: Gain with 
uni-directional pump 

Let us now consider the cascade/ladder scheme (as 
shown in Fig. |2|b)) with a uni-directional pump from 
the lower levello) to the upper level |c) at a rate k. To 
the zeroth order approximation in the probe field we 
obtained for the steady-state populations 



(0) _ KiYaTca + 2Q^) 



(0) _ YbiY'i^ca + 2Qg) 

&bb <Jl 



(Cl) 



(C2) 



(C3) 



where % = xiyb^ca + 2Q^J + {yb + K)(y„r„ + 2Q2). Here 
we also assumed resonant drive and probe excitation. 
In the strong field limit » yb,ya,K, we obtain the 
non-inversion condition as 



K<yb 



Also for gain we obtain 



K > y/YbYa 



(C4) 



(C5) 



Combining Eq.|C4 C5| , the condition for gain without 
population inversion gives 



(C6) 



iYaYb <K<Yb, and Yb > Y^ 



We see from Eq.| |C6[ that the spontaneous emission rate 
on the lasing and drive transition is opposite to that in 
Lambda-configuration. 
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l<3l> ■ 



\b)- 



FIG. 14: Three-level atom in which the transition \ai) —> \b} and 
\o2) —> \b) are coupled to the vacuum modes 



Appendix D: Vacuum-induced coherence 

In this section we will briefly review interaction of 
a three-level system with the vacuum. The three-level 
atom has two closely spaced upper levels as shown in 
Fig.|[l4|. The interaction Hamiltonian in the interaction 
picture is 



(Dl) 



Redefining the coupling constant as gk,a{^) = gk,a^'^^, 
we write Eq. ||D1| as 



H.c] 



(D2) 



where the coupling constant g'^ ^ and the atomic transi- 
tion operator 0,y are defined as 



,/ _ Pi) ■ ek,rt£k,a 



n 



(D3) 



(D4) 



Here again we will assume (without the loss of gener- 
ality) that the atom in the upper level at f = and 
the field modes are in the vacuum state |{0}). The state 
vector at f > is given by 

mt)) = c„,(f)k)l!0}> + c«,(0l«2>l(0})+ 

CiMAt)\mk,a) 

with Cfli(O) = l,Cfl, = and Ck,«(0) = 0. The evolu- 
tion equation for the state vector is governed by the 
Schrodinger equation iH\W{t)) = ^|^(f)) which gives 

C«, (t) = -i g|:;j'\R)e'<"">-^'-)'Q,k,„(f) (D6) 

k,a 



C«.(0 = -i X ^£f(R)e'^""^^"'^^'Qk,,(0 (D7) 



k,a 



-fg;<;f'(R)e-'K-n.).C„,(0 
Solving for Ci,^k,a(f) we obtain, 

Jo 

Jo 



(D8) 



(D9) 



Formal substitution of Eq.( |D9| in Eq.|D6 1 we obtain 

kj ' ' -'o 
Using simple algebra, we can obtain 



(DIO) 



• r e'l("»if.-n)f-Ki,-n-)f']Q^(r)rfr 
Jo 



(Dll) 



where the factor p is the most important dipole- 
orientation term defined as 



p = 



|p(flifc)||p(fl2fc)| 



(D12) 



Let us consider that the splitting between the two levels 
CO12 = (Oaib - oj^h ^ o^mbfC^'a^b and using the Weiger- 
Weiskopf and Markov approximation approximation] 2 1, 
we obtain 





; j g'[('U«i'-<i'a2t'')g-n't(t-f')j^f' 

Jo 



(D13) 



Performing the integration over the frequency variable 
and using the result e-'^'^^'-''^dvk = 5{t - t') we obtain 



Q(f) = -yG,(f)-p^^C«,(f) 

I e'[K'-"»2f'')6(f-f')dr 
Jo 



(D14) 



Here we have neglected the Lamb-shift term. Thus we 
arrive at the final form the Eq.| D14| as 



C.,(f) = -|c.,(i) - (D15) 
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FIG. 15: Two-level model. The decay rate from the levels \a) 
and \V) is giver\ by y while the rate of populating upper level 
is r. 



The second term is counterintuitive and results from the 
coupling of the two transition with the vacuum modes. 
Similarly we can obtam the evolution equation for Qjlf) 
as 



(D16) 



To eliminate the fast oscillating term e""'^^^ we will make a 
transformation Ca^{t) = ai{t) and Caj(f) = a2{t)e~'"^^' , thus 
Eq.( |DT5) ,l [DT6| l gives 



Hi) = 



(72 . \ 



ai{t) 



(D17) 



(D18) 



If we consider, degenerate upper levels then Eq.lDTSj 
gives 



Thus even vacuum of electromagnetic field can induce 
non-zero coherence between degenerate levels. 



Appendix E: Density Matrix Vs Rate Equations for TLA 

We consider a two level system with \a) and \b) as the 
upper and the lower levels. Population from the two 
levels decay out at a rate y while the rate of populating 
upper level is r. The density matrix equations are given 
as 

Qm = r- yQaa - 1 Qab - (El) 

Qbb = -yQbb + i {^'9ab - ^e'ai) (E2) 



.■a 6 
a 

D 

^4 



. 1 ' ' ' 1 ' 


1 1 ■ 1 1 1 ■ ■ 1 1 ■ ■ 1 1 ■ ■ ■ 1 . 

Inversion 




Photon 






■/ 





6 

t(ns) 



10 



12 



FIG. 16: Numerical simulation using the Rate equations. Using 
the parameters fc = 30 ns"^ . 



Qab - -TabPab - (.Pan - Qbb) 



(E3) 



The propagation equation for the field (Q) in the slowly 
varying amplitude approximation as 



where the coupling constant rj is 

rj = vNp^/2eocH 
Substituting p„;, = in Eq.dESl we obtain. 



Qab - -2 7;— {Qaa " Qbb) 
i-ab 



(E4) 
(E5) 
(E6) 



Substituting Eq.||E6| in Eq.ljETJ and Eq.|E2 1 we obtain, 

2Q2 



Qbb = -yQaa + 



(D19) obtain. 



- (Qaa - Qbb) 
(Qaa - Qbb) 

4Q2 



Qaa - Qbb - 1' ~ y (Qaa - Qbb) ^ — (Qaa - Qbb) 

i-ab 



(E7) 
(E8) 

(E9) 



Assuming the spatial uniformity of the field and using 
Eq.||E6|, the propagation equation for the field gives. 



-^ = ^{Qaa-Qbb)0. (ElO) 
at I ah 

Using simple algebra we obtain (for real field) 

dQ2 2a] , 
-rr = ^ (Qaa - Qbb) C2 (Ell) 
at Tab 

From the definition of Rabi frequency Q and field am- 
plitude £ we can write. 



Q = p&/2n, £2 = nhv/eoV 



(E12) 
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Here n is the number of photons. From Eq.lEll I and 
Eq.( E12| we obtain. 



dn p^vN 



{Qaa - Qbb) n 



From Eq.||E9| and Eq.|[ET2) we obtain. 



(E13) 



Qaa - Qbb - ^ - y {Qaa - Qbb) ~ 



eo^'^ab 



n (Qaa - Qbb) (E14) 



Let us define some parameters to get our result in con- 
sistent withl67J 



N„t = NV {q„„ - Qtb) , Rp = NVr, K = 



v9 



(E15) 



Using new parameters, our equations takes the form 

dn 

— = ¥Ma,n, N„b = Rp - yN„b - KnN„b (E16) 



To work on the numerical simulations we can use the 
rate equations derived from the density matrix equa- 
tions for different choices of T2. We have also added 
the cavity decay term in the equation of motion for n 
phenomenologically we gives us. 



dn 
It 



p^vN 
eoK^abV 



(Qaa - Qbb) n - Yen 



(El 7) 



In Fig 16 we have plotted the population inversion and 
number of photons generated. We can clearly see the 
relaxation oscillations. To summarize, rate equations 
can be derived from the density matrix or often known 
as master equation when the relaxation rate of the off- 
diagonal term Qai, is high thus reaching steady -state value 
much faster than any other time scale.. Similar equations 
can also be easily obtained for three-level and the four- 
level model for traditional lasing schemes. 
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